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Abstract 

We have used scanning gate microscopy to explore the local conductivity of a current-annealed 
graphene flake. Mapping the local neutrality point (NP) after annealing at low current density 
reveals micron-sized inhomogeneities which reflect the temperature distribution generated by ohmic 
heating. Broadening of the local e-h transition is also correlated with the spatial homogeneity of 
the NP. Annealing at higher current density improves the NP homogeneity, but we still observe 
some asymmetry in the e-h conduction. We attribute this to a hole doped region close to one of 
the metal contacts combined with underlying striations in the NP distribution. 



PACS numbers: 



Graphene exhibits a wealth of properties relevant to a wide range of applications and 
fundamental research.^ To meet expectations and complement silicon in future nanoelec- 
tronics it is necessary to obtain precise and consistent control over its nanoscale electronic 
properties. Promising steps in this direction have been made using a combination of chem- 
ical functionalization and geometrical confinement J 2 l 3 l 4 l S El Developing techniques for ana- 
lyzing and controlling the effect of these processes is thus at the forefront of graphene 
research. Scanning probes in particular have provided valuable insights into the substrate's 
influence on nanometer-scale topographic and potential fluctuations) 7 1 8 1 9 1 1 ^ At the mesoscale 
photocurrent! 11 1 12 1 and Ramari ^ 1 14 1 15 1 16 1 microscopy have also been used to characterize dop- 
ing from charged surface impurities and charge transfer from metal contacts. The presence 
of charged impurities manifests in the electrical tranport as a shift in the charge neutral- 
ity point (NP)p^ but the effect on the mobility and minimum conductivity at the NP is 
still debated Mesoscopic (micron-sized) inhomogeneities in the impurity density also 
contribute to the random variations observed in the transport properties of as-prepared two 
terminal devices. 20 Annealing in an inert atmosphere and degassing reduces the NP shift and 
inhomogeneitypEl b u t unless a final cleaning step is performed in situ the improvement is 
limited due to the re-adsorption of atmospheric gases and water vapor pZl Current-annealing 
has been used for in situ removal of adsorbatesj 23 1 24 1 in many cases resulting in an electron- 
hole conduction asymmetry. | 2 | 23 | 25! Here, we use scanning gate microscopy (SGM^to explore 
the local conductivity in a current- annealed graphene monolayer and show that persistent 
inhomogeneities in the remanent impurity density contribute to the anomolous e-h conduc- 
tion. We investigate a graphene flake (~ 8 x 8 fim 2 ) mechanically exfoliated from natural 
graphite onto a highly doped Si substrate capped with a 300 nm thick Si02 layer. Optical 
microscopy^ and Raman spectroscopy^ W ere used to locate the flake and confirm that it is 
a monolayer. Two 50 nm thick Au contacts were patterned using e-beam lithography and 
lift-off processing [see inset of Fig. [jja).] The sample was annealed at 200° C in N 2 /H 2 (5%) 
to remove resist residue and mounted to the head of a scanning probe microscope evacu- 
ated to 10 -5 mbar. Further cleaning was performed in situ by driving J ~5xl0 8 A/cm 2 
through the device and monitoring the shift in the neutrality point voltage (Vnp)- 23 The 
two terminal resistance R as a function of voltage Vbg applied to the Si back-gate [Fig. 
[lja)] reveals shifts of ^3 V in Vnp every 10 minutes. We assume the initial Vnp ~25 V is 
due to hole-doping adsorbates such as H 2 and 2 . Using n—a{VBG — Vnp) {ol—7.2 x 10 10 
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FIG. 1: (a) R(Vbg) measured at 10 min intervals (blue, solid) (arrow indicates the direction of 
shift in peak position,) and after 9 hours (red, dotted) of current annealing at J ^5xl0 8 A/cm 2 . 
After 10 mins at J ^lxl0 9 A/cm 2 (green, dashed). Inset: False-color atomic force micrograph of 
the sample (scale bar 2 jim). (b) Schematic of the setup used to perform SGM. 

cm -2 ) 17 to estimate the carrier density, we obtain a mobility of ~8000 cm 2 V _1 s _1 at n ~2x 
10 11 cm -2 . The final cleaning current was applied for ~9 hours, resulting in a shift in the 
bulk Vnp from ^25 V to ^16 V [red curve, Fig. [jja)], and a pronounced shoulder appearing 
at Vbg— 6 V. The latter is characteristic of flakes divided into regions with different carrier 
density, either by chemical dopingpES or invasive metal contacts.^ The suppression of R 
is also typical of flakes with a mesoscopic NP inhomogeneity. 20 Our SGM setup is shown 
in Fig. [jjb) (see Ref.[26] for details.) To benefit from the high signal-to-noise ratio achiev- 
able using a.c. detection, we modulate the potential difference Vt between the tip (Pt/Ir 
coated NanoWorld ARROW-NCPt) and the graphene at low frequency (typically 3 V @ 1 
kHz) and detect the modulation of Ids (~250 fiA) using a lock-in amplifier. We quantify 
the demodulated component (5Ids~0-1 mA) by the local transconductance g m = dlns/dVr 
normalized to the bulk conductance G(Vbg)^ 

Figure [2](a) plots g m (VBG) measured with the tip in each of the regions indicated in the 
diagram of the device (inset, Fig. [2|) G is suppressed (g m < 0) at each tip position when 
Vbg <6 V. The opposite is true for Vbg >16 V, when G is everywhere enhanced (g m > 0). 
The crossover back-gate voltage Vo = Vec^m— 0) falls within the range of intermediate 
Vbg and shows a sensitive dependence on tip position. To understand the basic form of 
the measured ^(Vbg)? we relate the conductivity (Jl of the ^(50 nm) 2 region perturbed 
by the tipP to a local carrier density by a^x^y) = \oVbg + PVr + n\e/i : where j3 is the 
capacitive coupling between the tip and graphene, and n is the local impurity-induced charge 
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FIG. 2: (a) Normalised transconductance measured with the tip at the locations marked in the 
inset, (b) Contour plot of the neutrality point Vq of the flake in the states corresponding to the 
red (1) and green (2) back-gate sweeps in Figjlja) (Scale bar 2 /mi.) 

density.^- 9 1 30 1 31 1 Neglecting quantum contributions to the conductivity, g m is proportional to 
da^x, y)/dVr, which changes sign at the local NP when Vbg — Vo ~ n(x,y)/a (since 
(Vt)=0). Hence Vb is a measure of n, which can be empirically related to the impurity 
density rii via n oc n\ (b =1.2-1.3).^ We explore n(x, y) by measuring g m {VBG) with the tip 
positioned over a grid of points with pitch ^300 nm and constructing the map of aVo(x, y) 
shown in panel 1 of Fig. [2](b). A band with n ~ 4.3 x 10 11 cm -2 (6 V) runs across the flake 
parallel to the contacts, asymetrically flanked by regions with higher density close to the 
contacts, where n ~ 11 x 10 11 cm -2 (16 V). This n profile can be explained by assuming 
that it reflects the temperature distribution generated by ohmic heating during current 
annealing. The remanent n distribution in this case suggests that energy is dissipated in 
the flake itself, while the contacts act as heat sinks. A similar temperature distribution was 
inferred from previous scanning prob^Sl anc [ Ramaii^l microscopy measurements, though 
here the asymmetry either side of the higher temperature region is more pronounced. The 
two values of Vo also coincide with the positions of the shoulder and the maximum of the bulk 
R(Vbg) [ re d curve, Fig. [jja)], confirming that these features originate from the mesoscopic 
inhomogeneity. Figure |3](a) shows plots of the local conductivity o- l (Vbg) reconstructed by 
numerically integrating g m (VBG) [Fig. [2|a)]. Within the framework of Ref. [30], the width of 
the minimum conductivity plateau AVbg at the NP is expected to increase with n. | 19 | 3Q | While 
we did not observe this direct relationship, Fig. [3](b) reveals a linear correlation between 

4 




FIG. 3: (a) Local conductivity at the same locations marked in Fig. [2^b) (reconstructed by 
integrating g m (VBG)-) Each curve has been vertically offset such that 0l(Vo)=O. Correlation plots 
between AVbg an d AVb, the change in Vq along linescans in the vicinity of the dashed line in Fig. 
[2^b). Dashed line is drawn to indicate the linear trend. 

AVbg and the spatial gradient of Vb along parallel lines in the vicinity of the dashed line in 
Fig. [2](b). (We extract AVbg for each curve by subtracting the back-gate voltages where 
&l ~ 0.2.) This is understandable as broadening is most pronounced when the tip gates 
both n- and p-type regions, in the same way that AVbg m homogeneous flakes is broadened 
by e-h puddles at low carrier density.^ 

To assess the inhomogeneity remaining when V/yp ~ 0, we annealed the flake at J ~1 x 10 9 
A/cm 2 for 10 minutes and repeated the measurement of n(x,y). The resulting n map is 
shown in panel 2 of Fig. [2](b) and the corresponding R(Vbg) sweep in Fig. [j](a) (green 
curve). In line with the previous analysis, the single peak in R{Vbg) coincides with the 
average value of Vo (~ 2V) (n = lx 10 11 cm -2 ), and the greater homogeneity is reflected by 
the absence of the shoulder. However, upon close inspection we observe a region close to the 
top contact with n = 7x 10 11 cm -2 . The doping in this region, which persisted even after 
annealing with J >lxl0 9 A/cm 2 , is responsible for the remaining asymmetry in R(Vbg), 
since opposing changes in R occur in regions with opposite carrier typeP^ Figure |I](a) shows 
a raw g m (x,y) (Vbg= 3 V) image of the flake after the second current anneal [c.f. green 
curve, Fig. [1(a)]. Such images provide complementary information to the maps of n(x, y) 
by allowing one to resolve finer n variations directly via changes in g m . Superimposed on 
the micron-sized inhomogeneity are pronounced striations in g m that span the length of the 
flake with lateral period of ^100 nm (see linescan, Fig. [4](b)) . These modulations could be 
caused by enhanced heating from resistive hotspots in the disordered pn junctions formed 
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FIG. 4: (a) SGM image (12 /imx5 nm) mapping g m at Vbg= 3 V [c.f. green curve, Fig. [ij a)], 
(b) Profile of g m along the line in image (a). 

close to the metal-graphene interface, 20 or electromigration of material from the contacts.^ 
By inspecting the amplitude of the modulation in g m) we estimate that these striations 
reflect an impurity density fluctuation An ~ 7.2 x 10 9 cm -2 , which may impose an intrinsic 
limit to the homogeneity achievable when current annealing supported flakes. 

In conclusion, SGM is a powerful method for characterizing the local conductivity of inho- 
mogeneously doped graphene. We find that the local NP reflects the temperature distribu- 
tion generated by ohmic heating and also exhibits finer linearly correlated inhomogeneities. 
Both types of inhomogeneity are likely to contribute to the e-h asymmetry observed in 
current-annealed flakes. 

This work was financially supported by the European GRAND project (ICT/FET). ACF 
acknowledges funding from the European Research Grant NANOPOTS and the Royal So- 
ciety. 
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